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Abstract 

Wc introduce a modified version of the two-dimensional Navicr-Stokcs equation, preserving 
energy and momentum of inertia, which is motivated by the occurrence of different dissipation 
time scales and related to the gradient flow structure of the 2-D Navier-Stokes equation. The 

^ , hope is to understand intermediate asymptotics. The analysis we present here is purely formal. 

Oj, A rigorous study of this equation will be done in a forthcoming paper. 

3 ' 1 Introduction 



The two-dimensional incompressible Euler equation in vorticity formulation reads 

{dt + u-V)uj{x,t)=0. (1.1) 

Here x G M?, t G M+ and u = u{x, t) G M? is the velocity field defined as: 



^: u = V^ij, ^ = -A-^w. (1.2) 

OO 

^^ . Explicitely, we have : 

1 1 X"' 

This equation is formally hamiltonian (we refer for example to [MTWR]). The Hamiltonian is the 
energy 

E{uj) = - -ipujdx. (1.3) 

By Noether Theorem, the center of mass M = J xuj (which is related to the invariance of E with 
respect to the group of translations) and the momentum of inertia (which is related to the invariance 
of E with respect to the group of rotations) I{u>) = 2 / 1-^ ~ M\'^Lodx are also conserved. Moreover, 
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due to the degeneracy of the Poisson bracket, one of the main feature of (jl.ip is the presence of an 
infinite number of conserved quantities, sometimes called Casimir: all the integrals of the form 

F^{uj) = / 4>{uj)dx 

are conserved. For this reason, the rigorous study of the large time behaviour of the solution of 
(jl.ip and hence the justification from (jl.ip of the presence of coherent structures observed in real 
and numerical experiments (see e.g. [MS]) remains widely open. In the following, we shall focus on 
non-negative solutions, we shall assume that a; is a probability distribution so that oj is non negative 
and J" w = 1 constantly in time. We also fix the reference frame in such a way that J xoo = for 
every time. 

One attempt to justify the appearance of these coherent structure is due to Onsager [O], see 
also [LP], [MJ]. The main idea is to replace the Euler equation (jl.ip by the system of N point 
vortices and to study the Statistical Mechanics of these point vortices. In a limit A^ -^ +00, the 
Gibbs measure associated to the point vortices concentrates to some special stationary solutions of 
the Euler equation (called mean field solutions), this was rigorously justified in [CLMPl], [CLMP2], 
[K] and [KJ]. These states are under the form: 

u = ^-^ (1.4) 

where 



is a normalization to have Jo; = 1. Recalling that lo = —Aip, we realize that eq.n ()1.4p is a 
nonlinear elliptic equation. The study of this equation in connection with the variational principles 
arising from statistical mechanics was performed in [CLMPl], [CLMP2] . Let us summarize the 
main results. We define the free-energy functional as 

Fib,a) (^) = S{lj) - hE{uj) - aI{oj) (1.5) 

for a given pair a < and 6 > 0. -^(fe,a) is defined on the space F of all the probability densities on 
M^ with finite entropy, energy and moment of inertia. We define the canonical variational principle 
as 

F(6,a) = infF(,,,)(w). (1.6) 

Next, for i? G M and / > let us introduce the set 

^{EJ) = V G ^Eiu) = E, I{co) = 1} (1.7) 

and consider the microcanonical variational principle 



S{E,I)= inf S{io). (1. 



^<^^{E,I) 



The above variational problems can be related to the solutions of the Mean-Field equation (jl.4p . 
Moreover, in the whole space M^, we have also the equivalence of the ensembles: 

Theorem 1 (CLMPl, CLMP2) For a < and < b < Svr; 

i) There exists a unique, radially symmetric minimizer lo = ^(h,a) ^ T o/ the problem (|1.6p 
which is the unique radially symmetric solution to eq.n (|1.4p . 

ii) When b ^ Sir, to converges (weakly) to a 5 at the origin. 

Hi) F{b, a) is a concave smooth function and 

OF dF 

-Q^ = H^{b,a)), -Q^ = -E{uJ(h,a)) 

iv) For i? E M and / > define 

S*{I, E) = sup(F(6, a) + bE + al) 

a,b 

and denote by b{I,E) and a{I,E) the unique maximizers. Then S{I,E) = S*{I,E) and hence S 
is a smooth convex function. 

v) The variational problem (jl.Sp has a unique minimizer 10(1, E) and 

^i^i E) = ^{b{I,E),a{I,E))- 

Note that when 6 < the theory is easier. Indeed the functional Fff,^a)i'^) is convex so the mini- 
mization problem is standard and eq.n (|1.4p has a unique (radial) solution [GL]. We also point out 
that the equivalence between (jl.6p and (jl.Sp which is established in v) is also useful to establish the 
existence of a minimizer for (jl.Sp . Indeed, in an unbounded domain, the existence of a minimizer 
for (jl.Sp seems difficult to establish directly because of the absence of higher moment control which 
would allow to pass to the limit in the constraint for the moment of inertia. Note also that eq.n 
(|1.4p has a natural statistical mechanical interpretation, its solutions being Gibbs states with a 
self-consistent interaction. Therefore —b is an inverse temperature. Hence 6 > implies negative 
temperature states, as predicted by Onsager [O] in terms of point vortex theory. 

A rigorous justification of the fact that the solutions of the mean field equation plays a special 
part in the large time behaviour of the Euler equation still seems an out of reach problem. An 
attempt towards the justification of the fact that the states ()1.4p play a special part in the 2D 
turbulence could come from the study of the intermediate behaviour of the Navier-Stokes equation. 
Indeed, for the Navier-Stokes equation, 

{dt + u-S/)u;{x,t) = i^Au;{x,t), (1.9) 

J LO and j xio are still conserved so that we can still consider non-negative solutions so that f lo = 1 
and J XLO = 0. Nevertheless, due to the dissipation term in the right hand side of eq.n (jl.9p . the 
asymptotic behavior of the solutions is trivial, namely lo{x, t) ^ pointwise and in the L^ sense for 



p > 1. Consequently, one can hope to observe the mean field solutions only as intermediate states. 
To formalize this idea, let us notice that the momentum of inertia / increases by a constant rate: 

i{uj) = 2v (1.10) 

consequently it can be considered as constant for times << Xjv. In a similar way, E and S are 
dissipated with rates 

E = -ufu;^ S = -uf^-^^. (1.11) 



Looking at eq.ns (jl.lip . one realizes that the energy could also evolve on a different and longer 
scale of times with respect to 5 (whenever the last term in (jl.lip dominates on the first one). This 
would suggest to consider, in the first approximation, E and / as constant, by looking at a master 
equation which modifies the Navier-Stokes equation leaving constant both energy and moment of 
inertia, but retaining all the other features of the Navier-Stokes dynamics. The derivation of such 
an equation, based on geometric arguments is the aim of the following section. 

2 Derivation of the model 

A natural and fruitful way to approach the problem is to invoke a recent characterization of the 
Navier-Stokes equation connected with the mass transport problem and the associated differential 
calculus introduced in [Qt| . We follow the excellent monographies [V], |AGS| for outlining the main 
ideas. 

Let M be the manifold of the probability measures on M^. One can formally give to A^ a 
structure of Riemannian manifold. For any p £ M we parametrize the tangent space as 

TpM = {p\p = -divu}. (2.1) 

Eq.n (|2.1|) expresses the tangent space to any point p of 7W as mass preserving velocity vectors p. 
Next, we define a Riemannian metric. On the tangent space TpAi, we define a scalar product by 

{pi,P2)w= p^^ui-U2dx, (2.2) 



being pi = —divui, i = 1,2. The gradient Vvf with respect to this Riemannian metric of a 
functional F : A^ ^ M, is defined as: 



f 5F 

{VwF,p)w = DF ■ p = - j — divu, 



where DF is the differential of the map F and 5F/5p the usual variational derivative. An explicit 

computation shows that 

SF 
VvKF = -div[pV— ]. (2.3) 



In a similar way, we can define on TpM a skew-symmetric operator Jw by 

JwP = — div iu ) , p = — div u (2-4) 



where w^ = {u2, —ui). Note that this yields in particular the expression 

s:rp 

Jw^wF = -d\v\pV^—-] . (2.5) 

op 

Gradient flows with respect to a functional F are the solutions to 

dtp = -VwF. (2.6) 

and are dissipative in the sense that 

d ^ f ,^6F,^ 

_f=_y.iv-p<„ 

whereas hamiltonian flows defined by 

dtp = Jw^wF 

are conservative since 

dt 
Since 5E/5uj = ijj, we get from (|2.5p that 

div wV V = Jw^wE 

and hence, the Euler equation can be interpreted as an Hamiltonian flow in this framework. 
Next, to interpret the dissipative part of the Navier-Stokes equation, we notice that 

Aw = div ( i^V— J = —VwS 

where S is the entropy functional. 

In conclusion the Navier-Stokes equation can be expressed in terms of a gradient and an anti- 
gradient flow : 

dtuj = -vVwS + Jw^wE. (2.7) 

According to the previous discussion we assume that the energy and the moment of inertia are 
varying much more slowly than the entropy functional. Therefore it may be useful to derive an 
effective equation according to the following prescription. We shall take the orthogonal projection 
(with respect to the scalar product (|2.2p ) of the vector field in the right hand side of eq.n ()2.7p on 
the manifold E = const and / = const, with the aim to characterize the states which are close to 
the true dynamics for a large interval of times (coherent structures), as asymptotic states of the 



new dynamics. We shall see that such states are the solutions to the mean-field equation (jl.4p . 
Since: 

VwE{uj) = -div[u;V— ^] = -div[a;VV'], 

Vh//(w) = -divT^V— r^l = -divLV— 1 = -divlujx]. 
'- ouj '- 2 -' '- -' 

we find that 

(Vi^-E, Jw'VivE)]y = (VvK-^ ■ Jw^wE)w = Oi 

and hence Jw'^wE is tangent to the manifold E =const, / =const. On the other hand the 
projection of VwS on the tangent space of such a manifold is of the form 

ViyS — b div(ix'V^) — a div(ti;x), 

with a and b two suitable multiplicators. As a consequence, the equation we are looking for is: 

dtu; + u ■ Vlo = V div(Va; — bcoVip — aujx) (2-8) 

= V div\ujV {log oj — bijj — a— )\, 



with a and b to be determined by the simultaneous conservation of E and /. A straightforward 
computation yields: 



2//a;2 + 2y _ 2/cj|V^|2 + y/ 



uo 



2 



^ 2//cj|V^|2-y2' " 2lJuj\V^\^-V^ ' ^^'^^ 

where 

V = u!x ■ Vip = dx dyuj{x)u!{y)x ■ Vg{x — y) = — — ■ (2.10) 

Note that we have by the Cauchy-Schwarz inequality 

V^ = {[ ujx-Vi^f <2I I uj\Vil}\^, (2.11) 

and hence b is positive if inequality ()2.1ip holds strictly and J lo^ > ^^. 

We point out that eq.n (12. 8|) has been derived in [Ch2] (see also [Chi]), by using different 
argoments than thse of the present paper. 

Note also that another class of similar equations was introduced in [RS] and studied mathemat- 
ically in [MR]. Nevertheless, despite to some formal analogy the mathematical properties of the 
equations studied in [MR] are very different from the ones presented here. 

Let us now discuss what we may expect about the asymptotic behavior of eq.n (j2.8p . The main 
feature of eq.n (j2.8p is the decay of the entropy functional. Indeed, we have 

dS{uj) _ dS{uj) dE{uj) dl{uj) _ 
dt dt dt dt 



,SS JE Sr, r ^,5S 61 SE^ 
b- a^)div\io\/[- a- o-^j 

'OLO du> OUJ OLO du> du> 

-./u;|V(--a--6-)| = 



-ly / uj\V([ogLLi — bip — a 



I |2 



2 ' 

In particular, this suggests formally that the asymptotic states satisfy 

ujV{loguj — bijj — a ) = 

and hence the mean field equation p.4p . 

We remark that the procedure of constructing dissipative equations leaving invariant a given 
quantity is not unique. For instance the heat equation in M?: 

dfUJ = Alj 

can be modified to leave invariant / according to the procedure suggested by the gradient flow 
structure. The result is 

dtU! = div(Vw — auox), (2.12) 

where 

1 
a = -j. 

On the other hand we could also have 

dtio = dleuj. (2.13) 

Note that eq.ns (I2.12|) and (|2.13p have different asymptotic states. 

An attempt to characterize an intermediate asymptotics was presented by Gallay and Wayne 
[GWj according to the following ideas. 

It is well known that a special solution to eq.n (jl.Op (for u = 1) is given by the so called Oseen 
vortex: ^ 

Note that this is also a solution to the heat equation. It was shown in |GWj that this solution 

describes the long time asymptotic of the Navier-Stokes equation in L^. Indeed, with the change 

of variables 

^ = -^; r = log(l+t), uj{x,t) = {l + ty'w{^,T), 

V J- + * 

the Navier-Stokes equation in the new variables is under the form : 

drW + v-V(^w = A^w+ V^- (-^w). (2.14) 



It is possible to show that w — > M^ in L^ as t — > c«, where W{S,) is the rescaled Oseen vortex. As a 
consequence the Oseen vortex can be thought as characterizing an intermediate asymptotics before 
the dissipation scale. Note that W is also a solution to (|1.4p for b = 0. 

This analysis enters perfectly in the context of the projected gradient flows. Indeed neglecting 
the energy, the mere Constance of / yields 

dtu; + u-Vuj = Auj + -V ■ (ujx)), (2.15) 

that is eq.n (j2.14p for 1 = 2. In this particular case we have seen how eq.n (|2.14p can be obtained also 
by a simple change of variables, due to the fact that the dissipation rate of I is constant. Roughly 
speaking the analysis we present here is an attempt to see what happens before the appearence of 
the Oseen vortex by projecting on a less robust manifold. Indeed one could argue that / is more 
stable than E in many interesting physical situations. If so eq.n (j2.8p should be more appropriate 
on the time scale when E is practically constant, while eq.n (j2.14p should describe the fluid when 
E start to be dissipated at constant /. After that, everything disappears. 

It would be interesting to look for a numerical evidence of this fact, if true. 

The formal derivation of eq.n (j2.8p is also in agreement with the stochastic vortex theory as we 
are going to illustrate. 

As explained previously, the procedure of constructing dissipative equation leaving invariant a 
given quantity is not unique. Equation (j2.8p has been interpreted as a constrained Navier-Stokes 
flow and it turns out that the asymptotic states are solution of the microcanonical variational 
principle. Since this variational principle is obtained by Mean Field limit from the statistical 
theory of the point vortices, it is interesting to see how equation (|2.8p is related to the theory 
of point vortices. Now we show how, at least formally, the structure of the constrained gradient 
flow is compatible with the discretization of the Navier-Stokes equation obtained by mean of the 
stochastic vortex theory. Since in this context we are not interested in the asymptotic behavior, 
we limit ourselves in considering the energy constraint only. Also the viscosity does not play any 
role so we set v = \. 

Consider A^ stochastic vortices in M^. They obey the stochastic differential eq.n: 

dxi = — 2^ V g{xi — Xj)dt + V2dwi (2.16) 

3 

where {wi}f^i are N independent standard Brownian motions. Here g{x) is a regularization of the 
Green fuction — 2;^^ log |a;|. It is well known ( [MP2| . |0s] ...) that the empirical random measure 

1 ^ 
fij\f{dx,t) = —y^d{xj{t) — x)dx (2-17) 

i=i 
approaches the solution to the (regularized) Navier-Stokes equation, if it happens at time zero. 



We now consider the mean-field energy: 

H{xi ...xn) = j^^aixj -Xr). 

j<r 

In order to guarantee the condition 

dH = 

a short computation using the Ito formula shows that we have to modify the process according to: 

dxi = T7 X^[V"^5(^i - Xj) - bN{t)Vg{xi - Xj)]dt + V2dwi+ (2.18) 

j 



where 



|VF|4 \VH 



AH f diiN^g * iJLN 

ON\t) ~ ~ 



|VF|2 /d^^|V<7*^Ar|2' 

Note that the above expression makes sense only if we regularize g and this explains why we did it. 
We observe that an analysis on the size of the last two terms in (j2.18p shows that they should 
be negligible in the limit N — > oo. Thus we introduce the essential process defined by 

'^yi = J^Y \^^9{yi - Vj) - bN{t)Vg{yi - Vj)] dt + V2dwi. 
j 

Now if we replace fi]\[{dx,t) by uj{x,t)dx in the limit N — > oo, each process Xj approaches the 
nonlinear (in the McKean sense, see |Mc] . |MP2] . |0s| ) process solution of 

dy = [V^g * uj{y) - b{t)Vg * oj{y)] dt + V2dw. (2.19) 

where 

, , , f dx LoAq * uj 



jl dxio\SJg * Lo\'^ 

From eq.n (I2.19p . we also derive the backward Kolmogorov equation for the probability distribution 
LO of the proces y: 

dtu; + V ■{uio) =V ■{Vuj-bLoVg*uj). (2.20) 

where u = V g * uj. It is remarkable that such an equation is the microcanonical equation ()2.8p . 
for a = 0, if we replace g by the true Green function in eq.n (j2.20p . 
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